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> ■ ABSTRACT 

Tij" \ Poststarburst (also known as K+A) galaxies exhibit spectroscopic signatures indicat- 

■ ing that their star formation was recently quenched; they are candidates for galaxies in 
' transition from a star-forming phase to a passively-evolving phase. We have spectro- 

scopically identified large samples of poststarburst galaxies both in the Sloan Digital 
in ■ Sky Survey (SDSS) at z - 0.1 and in the DEEP2 Galaxy Redshift Survey at z - 0.8, 

' using a uniform and robust selection method based on a cut in H/3 line emission rather 

00 . than the more problematic [O II] A3727. Based on measurements of the overdensity of 

■ galaxies around each object, we find that poststarburst galaxies brighter than QAL*g 

at low redshift have a similar, statistically-indistinguishable environment distribution 
as blue galaxies, preferring underdense environments, but dramatically different from 
that of red galaxies. However, at higher-z, the environment distribution of poststar- 
burst galaxies is more similar to red galaxies than to blue galaxies. We conclude that 
the quenching of star formation and the build-up of the red sequence through the 
K-l-A phase is happening in relatively overdense environments at z 1 but in rela- 
tively underdense environments at z 0. Although the relative environments where 
quenching occurs are decreasing with time, the corresponding absolute environment 
may have stayed the same along with the quenching mechanisms, because the mean 
absolute environments of all galaxies has to grow with time. In addition, we do not find 
any significant dependence on luminosity in the environment distribution of K-|-As. 
The existence of a large K+A population in the field at both redshifts indicates that 
cluster-specific mechanisms cannot be the dominant route by which these galaxies are 
formed. 

Our work also demonstrates that studying poststarburst-environment relations by 
measuring the K-l-A fraction in different environments, as is common practice, is highly 
nonrobust; modest changes in the comparison population used to define the fraction 
can drastically alter conclusions. Statistical comparisons of the overall environment 
distributions of different populations are much better behaved. 

Key words: galaxies: evolution — galaxies: formation — galaxies: starburst — 
galaxies: high-redshift — galaxies: stellar content — galaxies: statistics 
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1 INTRODUCTION 

How bulge-dominated (early-type) galaxies form is a long- 
standing puzzle in galaxy formation and evolution theories. 
These galaxies have old stellar populations and no ongoing 
star formation. However, they did not all form at high red- 
shift: many emerged at relatively recent epochs, as evidenced 
by the double-to-quadruplc increase in their comoving num- 
ber density since 2 1 ( BcU ct al. 20oi lWillmer et al.ll2006l : 
iFaber et all l2007l : iBrown et all 120071 ') . It is poorly under- 
stood what processes shut off the star formation in their pro- 
genitors and possibly transformed their morphologies and 
dynamical structures. 

Poststarburst gal axies, first identified by 
iDressler fc GunnI l| 19831 ). have been hypothesized to 
be the direct progenitors of early-type galaxies, since their 
spectra show signs of recent quenching of star formation: 
they possess large populations of young stars but no 
ongoing star formation. These galaxies are also called K+A 
or E-l-A galaxies, as their spectra can be roughly decom- 
posed into a combination of a K giant star (or early-type 
galaxy) spectrum an d an A star spectrum (see Fig. [1]). The 



visual morpholo gy (|Blake et all 12004 iBalogh et alTl2005l: 
Yang et alll2008l ). Sersic index distribu tion I Quintero et al.l 



20041 ). kine matics (iNorton et all [1^), and the elemental 



abundances (|Gotdl2007l ) of poststarbursts also suggest that 
they are likely progenitors of early-type galaxies. Studies of 
poststarbursts aim to reveal how star formation is quenched 
and shed light on the formation of bulge-dominated 
galaxies. 

An important hint for solving these puzzles comes from 
galaxy environment studies. The environment distribution 
of poststarbursts can differentiate among possible quench- 
ing mechanisms, as in some models the quenching is caused, 
or partly induced, by the external environment of a galaxy 
(e.g., by mechanisms that are strong only in groups or clus- 
ters of galaxies). Comparing the environment distribution of 
poststarbursts to their potential progenitors, blue galaxies, 
and possible descendants, red galaxies, will reveal what the 
preferred environment for quenching is. 

Previous studies of the environment distribution 
of poststarbu rs ts ha ve found contra dic tory results. 
Dressier et ahl l[l999l), iPoggianti et al.l (|l999l '). and 



Tran et al. ( 2003 , l2004h used samples around interme- 
diate redshift clusters (mostly 0.3 < z < 0.6) and found a 
higher poststarburs t fraction in clusters than in the field; 
iBalogh et all l| 19991 ) did a similar study around X-ray lu- 
minous clusters from CNOCl at 0.18 < z < 0.55 but found 
the p oststarburst frac ti on in clusters is similar to the field ; 
while IZabludoff et"aLl (Il996l) ( LCRS ) . iBlake et al.l (|2004|) 



f2dF) iQuintero et al.l ll 20041) . iGotd l|2005l ). IBalogh et all 
l|2005l) , and lHogg et"all (|2006l ) (SDSS in the last four) used 
local samples (all with z ~ 0.1) and found poststarbursts 
are more likely found in the field, outside clusters. These 
past studies, except those using 2dF or SDSS, usually 
were limited in sample size, especially at higher redshifts. 
The methods for selecting K-l-A galaxies in these papers 
also differed both in the spectral features adopted and the 
threshold used. The poststarburst sample s in m ost of these 
studie s, except for lOuintero et al.l (|2004l ) and iHogg et aP 
(|2006l ). are often incomplete or heavily contaminated due 
to selection on [O II] equivalent width, causing the large 



population of poststar b urst AGN hosts to be omitted, as 
explained in I Yan et al.l (|2006l ) and in Sec. IS.ll of this paper. 
Therefore, the poststarburst environment distribution is 
poorly known at high redshift [z > 0.2). However, high 
redshift is particularly important in understanding the 
formation of poststarburst galaxies since their fraction 
among all galaxies is significantly higher at z ~ 1 than 
2 ~ (Lc Bo rgne et al.ll2006l : I Wild et al.l | 2009l ). Whether 
the environments of poststarbursts evolve with redshift is 
also unknown. 

Most of these previous studies compared the fraction 
of poststarbursts among cluster galaxies with the fraction 
among a field sample. The result from this method depends 
sensitively on the parent sample in which the fraction is mea- 
sured. Because galaxy properties differ significantly across 
different environments, the difference between poststarburst 
fraction across different environments can be largely due to 
the different mix of galaxies in the sample. A meaningful 
comparison is to count poststarburst galaxies relative to 
their potential progenitors. A varying fraction across differ- 
ent environments for poststarburst relative to their potential 
progenitors would reveal how environment affects the possi- 
bility for quenching. 

In this work, we select poststarbursts uniformly from 
two large redshift surveys: SDSS at z ~ 0.1 and DEEP2 
at z > 0.7, using improved selection methods. The sample 
from DEEP2 is the largest clean poststarburst sample to 
date at high z. With these samples, we investigate the en- 
vironment distribution of poststarbursts at both low z and 
high 2, and compare with their potential progenitors and 
possible descendants. 

This paper is organized as follows. In Section 2, we de- 
scribe the data used. In Section 3, we describe the selection 
of poststarburst galaxies, the parent sample construction, 
the fraction and contamination estimates, and the environ- 
ment measurements. We then present the colour-magnitude 
distribution of poststarburst galaxies in Section 4. In Section 
5, we show the environment distributions of K-l-A galaxies 
in SDSS and DEEP2, compared with blue and red galaxies. 
We compare our results with previous studies and discuss 
the evolution with redshift. In Section 6, we investigate the 
luminosity dependence in the K-l-A environment distribu- 
tion. We discuss the redshift evolution of K-f-A environment 
distribution in Section 7 and the implications of these re- 
sults for different quenching mechanisms in Section 8. We 
conclude in Section 9. 

Throughout this paper, we use a ACDM cosmology 
with Q.M, Oa=(0.3,0.7), and a Hubble constant of H = 
lOOh km s-i Mpc~ ^ All mag nitudes within this paper are 
on the AB system (|Oke fc Gunn„198a ). 



2 DATA 



The DEEP2 Galaxy Redshift Survey (|Davis et all |2003| . 
|2007| . Faber et al. in prep) was a three year program using 
the DEIMOS spectrograph on the Keck II telescope cover- 
ing ~ 3 sq. deg of sky in four widely-separated fields, to 
a limiting magnitude of Rab = 24.1. The photometry was 
obtained with the CF H12K camera on the Canada-France- 
Hawaii Telescope (see lCoil et al.|[2003 ). To focus the spec- 
troscopic observations on high-z galaxies, in three of the 
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Figure 1. Sample poststarburst galaxy spectrum taken from the DEEP2 galaxy redshlft survey. The black curve shows the observed 
spectrum with a smoothing window of 11 pixels. The red and blue curves show the spectral decomposition. The templates used are the 
same as in lYan et al.l \200Qi\ . The blue curve is the young stellar population component, which is dominated by A stars. The red curve is 
the old stellar population component, which is dominated by K— giant stars. The green curve is the linear combination of the red and blue 
curves and it follows the data closely. This demonstrates that the K+A model fits the spectrum fairly well. Note that the [O III] A5007 
emission is prominent in this spectrum. The high [O III]/H/3 ratio indicates this poststarburst galaxy is also a narrow-line Seyfert galaxy. 



four fields, a colour selection is applied to select galaxies 
at 2 > 0.75. Tests done on the fourth field prove this pre- 
selection is successful: only ~ 10% of the final sample is 
from z < 0.75, while < 3% of actual objects at z > 0.75 
are discarded. Each CFHT pointing is adaptively tiled with 
40 slitmasks in an overlapping chevron pattern, giving each 
galaxy two chances to be on a mask. Tests done with mock 
catalogs i ndicate that dense regions are only slightly under- 
sampled l|Gerke et al.ll2005l ). Overall 60% of galaxies that 
meet the selection criteria are targeted. 

We have obtained spectra for ~ 50, 000 galaxies, with 
a resolution of i? ~ 5000 covering the wavelength range 
from 6500A to 9100A. The data are reduced by a dedi- 
cated reduction pipeline (Cooper et al. in prep). Redshifts 
are measured by the pipeline and confirmed by eye; > 33, 000 
galaxies have yielded high-confidence redshifts (> 95% con- 
fidence). Absolute B magnitudes and restframe U — B colour 
are computed using the K-correction code described in 
IWillmer et all (|2006l i 

The SDSS (|York et aklboOOl : IStoughton et al.ll2002l ) is 

an imaging and spectroscopic survey that covered n 
steradian of the celestial sphere, utilizing a dedicated 2.5- 
m telescope at Apache Point Observatory. The imaging 
are collected w ith five broadband fi l ters in drift scan mod e 
{u, g,r,i, and z: lFukugita et al .1 19961 : IStoughton et al.ll2002h . 
Spectra are obtained with two fiber-fed spectrographs, cov- 
ering the wavelength range of 3800-9200A with a resolution 
of R ~ 2000. The SDSS fibers have a fixed aperture of 3". 



The SDSS spectroscopic data used here have been 
reduced through the Princeton spectroscopic reduction 
pipeline (Schlegel et al. in prep), which produces the fiux- 
and wavelength-calibrated spectraQ The redshift catalog of 
galaxies used is from the NYU Value Added Galaxy Catalog 
(DR4) llBlanton et al.ll2005l]. K-c orrections for SDSS were 
derived using Blanton et al.l (120031 ) 's kcorrect code v3_2. 



3 METHODS 



3.1 Poststarburst Selection 



In principle, poststarburst galaxies are defined to have no 
ongoing star formation but to have experienced a significant 
star- forming epoch in the recent pastlj Practically, the two 
criteria have to be measured from specific spectral signa- 
tures, which unavoidably brings along the potential for bias, 
contamination, and a mixture of different physical processes. 

The first criterion, the lack of ongoing star formation, 
is usually indicated by the lack of line emission. However, 
line emission can also result from AGN (Active Galactic 



^ http://spectro.princeton.edu/ 

^ http:/ /wassup. physics. nyu.edu/vagc/ 

^ An important caveat regarding the necessity of a burst will be 
discussed further at the end of this section 
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Nuclei) or LINER(Low-Ionization Nuclear Emission-line Re- 
gionj3 activity. Therefore, a selection based on the lack of 
line emission will not be complete, as it will be biased against 
poststarburst galaxies with strong AGN/LINER-like emis- 
sion. The level of incompleteness caused by such a selec- 
tion method highly depends on the choice of emission line 
used. The two commonly used star-formation indicators in 
the optical are [O II] A3727 and Ha. Most previous K-l-A 
samples employed [O II], which is more easily accessible 
observationally at z > 0.4. However, as demonstrated by 
lYan et al.1 ()200 6l). significant [O II] emission associated with 
an AGN /LINER is observed in nearly half of the red galaxies 
and > 70% of K-l-A galaxies (as defined using Ha) , causing 
a large incompleteness in K-l-A selection. In contrast. Ha 
emission in LINERs is much weaker; they exhibit a very 
high [O II]/Ha ratio. In addition, [O II] emission from star- 
forming galaxies is much weaker than Ha due to dust ex- 
tinction. The combined effect makes poststarburst galaxies 
totally indistinguishable from star-forming galaxies accord- 
ing to [O II] emission, but clearly distinguishable according 
to Hq emission. Therefore, K-l-A selection is much cleaner 
if Ha is used as the star formation indicator. At z > 0.4, we 
have a dilemma as [O II] is so undesirable, while Ha is red- 
shifted into the infrared. We therefore rely on H/3 emission 
as a star formation indicator in this work; it is observable 
in optical spectra for z < 0.9. Although it is much weaker 
than Ha, it still provides a much mo re complete and more 
robust selection of K-|-As than [O II] l|Yan et al.ll2006l '). 

The second criterion for a galaxy to be a poststarburst, 
significant recent star formation, can be indicated by the 
presence of a large population of A stars, due to their short 
lifetime (< IGyr). These stars will cause strong Balmer ab- 
sorption features to be detectable in a galaxy's spectrum. 
Most past stud ies used the equival e nt width (EW) of H^ 



in absorption (iDressler et al 



iBalogh et all 1 19991 : iGoto et al.l |2003| ) , or a combination of 
two or three Balmer absorption lines ijZabludoff et al.l 1 19961 : 



19991: iPoggianti et al.1 [l999l : 



iTran et al.ll2003l . |2004| : iBlake et al.ll2004l ), to indicate the rel- 
ative l evel of recent st ar formation. Recentlv.^Q.uint ero et al.l 
l|2004h and I Yan et al.l (j2006) decomposed each galaxy spec- 
trum into two components - a young component mimicking 
A stars and an old component similar to K stars - and used 
their ratio (dubbed A/K) to indicate the relative strength 
of the recent burst. We take the latter approach. However, 
for robustness, instead of taking the ratio between the am- 
plitudes of the two components, we use the flux ratio be- 
tween the young component and the total stellar continuum 
around 4500A, denoted by /a for "A-star fraction," to in- 
dicate the relative level of recent star formation. This has 
several advantages over both EW indicators and the A/K 
ratio: it effectively averages over multiple Balmer absorp- 
tion features, providing a less noisy measurement than the 
average of two or three EWs; it is more stable than the A/K 
ratio as the latter becomes noisy when the K component is 
weak; and it is less sensitive to the overall continuum shape 



* It is still uncertain whether LINERs are associated with 
AGN. There are other non-AGN mechanisms that could produce 
LINER-like emi s sion li ne ra tios, none of whic h is tied to star for- 
mation. See ll2004h and lYan et al.l ||2006D for further discus- 
sion. 



or flux calibration as a broadly-smoothed continuum is re- 
moved before fitt ing. For details of the fitting procedure, see 
lYan et al.1 l|2006D . 

We therefore measure the two indicators, H/3 EW (in 
emission) and /a, in each galaxy spectrum from DEEP2 and 
SDSS. K-l-A galaxies are identifiable in a plot of H/3 EV>^ 
vs. /a. As shown in Fig. (2] most galaxies exhibit a tight cor- 
relation between these parameters. This is to be expected, 
as galaxies with relatively strong recent star formation gen- 
erally have a high ongoing star formation rate. K-l-A galax- 
ies stand out as a horizontal spur in this diagram, as they 
have no ongoing star formation but significant recent star 
formation. This separation between K-l-A ga laxies and star- 
formin g galaxies was first demonstrated by iQuintero et al.l 
(|2004h . Here we use an improved recent star formation indi- 
cator, /a, in place of the A/K ratio used in lQuintero et ahl 
(2004), for reasons mentioned above. This makes the sepa- 
ration much clearer, even though H/3 is weaker and harder 
to measure than Ha. The solid lines in the plots indicate our 
selection criteria for K-l-A galaxies, which can be expressed 
as the two inequalities: 



/a > 0.35 

and 

H/3 EW(A) < 4 X /a - 1. 



(1) 
(2) 



The choice of the /a threshold has its arbitrariness since 
the horizontal spur in Fig. [2] is a continuous sequence from 
strong K-l-A galaxies extending to quiescent, old galaxies at 
(0,0). Ideally, we would put the threshold at the point above 
which the K-l-A sequence (the horizontal spur) is separable 
from the star-forming sequence. However, this depends on 
the intrinsic scatter and measurement error of each popu- 
lation. In the SDSS panel of Fig. [2] the separation can be 
achieved at least down to /a of 0.25. But in the DEEP2 
sample, due to larger measurement errors, we are forced to 
use a higher threshold to keep the sample clean. Therefore, 
we choose /a > 0.35 as the criterion for both DEEP2 and 
SDSS samples in order to make fair comparisons. The H/3 
cut is chosen to go through the valley between the two se- 
quences in SDSS. Our results are less sensitive to variation 
in this cut as the valley is underpopulated. 

An important caveat in the definition is that even strong 
Balmer absorption features, such as H(5 EW(inabsorption) < 
— 5A, or equivalently /a > 0.4, do not guarantee the ne- 
cessity of a starburst, or a significantly elevated SFR, be- 
fore quenching. A sharp quenching of star formation with- 
out ab]n^t_is_^soj20^ihle_Jp produce a K-l-A galaxy. Us- 
ing [BruzuaT&^^harlo^ (|2003l ) population synthesis models, 
we found that a 5 x IO'^'^Mq (stellar mass) galaxy with 
a lOMo/yr constant SFR would have an /a of 0.6 (H5 
EW~ — 7A) after abrupt quenching, which does not drop 
below 0.3 for another 0.5Gyr. Such a level o f star forma- 
tion is common at redshift 2 ~ 0.7 and above (|Noeske et al.l 
I2OO7II . Therefore, a starburst is not absolutely necessary for 
a galaxy to be classi fied as a K-f A galaxy. This is also shown 
bv lLe Borgne et al.l (,2006.) using PEGASE-H R stellar popu- 
lation synthesis code ( Le Borgne et al1l2004h . Although the 
optical and NIR colours of K-fA galaxies suggest enhanced 

^ Throughout this paper, we use positive values for emission line 
EW and negative values for absorption line EW. 
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recent star formation is necessary (jBalogh et al.l 120051 ). the 
accuracy of the name 'poststarburst' is questionable. The 
name "K+A" is still precise since it only refers to the obser- 
vational definition, but it lacks a reference to the physical 
processes. A more appropriate name for these galaxies would 
be "post-quenching galaxies" . However, for consistency with 
the literature, we have been refering to all K+A galaxies as 
"poststarburst" galaxies. From now on, we will adopt the 
name "K+A galaxy" and "post-quenching galaxy" . 

3.2 Parent Sample Definition 

3.2.1 DEEP2 

The DEEP2 sample we use is a subset of the spectra ob- 
tained by the survey. We require galaxies to have H/3 cov- 
ered in the spectrum and to not be severely impaired by sky- 
subtraction residuals. The spectra are also required to reach 
a minimum restframe wavelength below 3900A to make re- 
liable spectral decomposition possible. This limits us to the 
redshift range 0.65 < z < 0.88, yielding a sample of 9564 
galaxies. However, we cannot use them all for the following 
environment analysis, for a few reasons. 

First, the DEEP2 sampling density drops precipitously 
below 2 ~ 0.7 due to the colour preselection of the survey, 
leading to undesirably large correction factors and high un- 
certainty in environment measurements, we limit the sample 
to z > 0.72. This removes 14% of the sample. 

Secondly, the environment measurement is not reli- 
able when a galaxy is too close to a survey edge. As a 
result, we further limit the sample to only galaxies that 
are at least I fe'^^Mpc (comoving) away from a survey edge 
(see ^3.41 and ICooper et al.l |2006| ) . This criterion removes 
23% of the sample after the above cut. In total, we have 
(1 - 0.14) X (1 - 0.23) = 66% of the original sample left. 

Thirdly, we select a luminosity-limited sample for 
the analysis. DEEP2's R-band magnitude limit {Rab < 
24.1) corresponds to a col or- and redshift-depe ndent limit 
in restframe B-band (see IWillmer" et aLlbOOeD . This cut- 
off is a function of color, shallower for redder galaxies 
l|Willmer et al.ll2006l ). At z = 0.88, the high redshift end 
of the K+A window, the limit is Mb — 51og/i ~ —19.7 for 
typical red galaxies with restframe U — B = 1.2. There- 
fore, DEEP2 is efi'ectively volume-limited for galaxies of 
all colours to Mb - 5 log /i < -19.7 for z < 0.88. With 
this magnitude cut, 53% of the sample after above cuts 
is further removed, most of them relatively blue galaxies. 
In total, we have 31% of the original sample left. We also 
make a bright subsample using a brighter luminosity cut at 
Mb — 5 log h < —20.7, which is approximately L* at z = 0.8, 
to investigate the possibility of luminosity dependence. 

In principle, to make an unbiased comparison between 
post- quenching galaxies and their potential progenitors and 
descendants, it is most ideal to make a stellar-mass lim- 
ited sample. However, we do not have the luxury of do- 
ing this. Altho ugh K-hand data is available in 1/3 of the 
DEEP2 fields l|Bundv et al.l 120061 '). the combined depth of 
the DEEP2 spectroscopy and the K-hand survey is not 
deep enough to cover enough post-quenching galaxies to al- 
low significant detection of any environment trend. There- 
fore, we resort to a B-band luminosity-limited sample. Our 
Mb — 5 log ft < —19.7 limit corresponds approximately to a 



stellar mass limit of M* > 10 ' MQ for red galaxies and 
M, > 10^°-^M© for blue galaxies. For K+A galaxies, it is 
approximately M* > 1O^"'®M0. 

Lastly, we limit our sample to only galaxies with a high 
signal-to- noise (S/N) spectrum. When the S/N per pixel in 
the continuum falls below 1 (corresponding to S/N~ 2.3 
per FWHM), the spectral decomposition becomes unstable, 
which also introduces large errors on the H/3 EW. This will 
contaminate the K+A sample with blue, star-forming galax- 
ies. We therefore remove all galaxies with S/N per pixel less 
than 1.0. This criterion removes 10% of the sample, most of 
which are blue galaxies. In addition, we remove a handful 
(0.5%) of galaxies which have H/3 EW error greater than 5A 
or have H/9 EW less than -5A(where we define emission as 
positive EW). These criteria effectively clean the K+A sam- 
ple without changing the whole sample significantly. They 
only brighten the median Mb by 0.06 for the whole sam- 
ple, 0.02 for red galaxies, and 0.06 for blue galaxies, all of 
which are insignificant. In total, we have 28% of the original 
sample left. Table [1] summarizes the details of our sample 
definition. 

With the above selection, we are left with 2649 galaxies 
from the DEEP2 survey. This is the primary DEEP2 sample 
used in the following analysis. There are 39 K+A galaxies 
in this sample. After inspecting their spectra individually, 
we manually removed 5 of them which have H/3 severely un- 
derestimated due to strong residuals from sky subtraction. 
The remaining 34 K+A galaxies are all robustly detected, 
brighter than our luminosity cut, and have well-measured 
environment. Whenever possible, we lift some of the selec- 
tion cuts and use a larger sample. For example, in Fig.|3] we 
lift the redshift cuts, edge distance cut, and the magnitude 
cut, but keep the signal-to-noise cut, giving the maximum 
sample of 74 K+A galaxies. All of the 74 K+As are plotted 
in Fig. [3] This is the largest K+A sample ever found in this 
redshift range. 

3.2.2 SDSS 

To test the evolution in the K+A population from z ~ 1 
to 2 ~ 0, we build a comparison dataset from the SDSS 
main galaxy sample in the redshift range 0.07 < z < 0.12. 
To make a fair comparison, we select a sample down to the 
same B-band luminosity limit relative to L», the character- 
istic luminosity of galaxies in a Schechter- function fit, as at 
z ~ 0.8. As star formation rates decline and stellar popu- 
lations age, galaxie s fade; L« in t he B-band decreases by 
1.3 mag per unit z (|WiUmer et al.! 2006; Fabcr ot al. 2007). 
Our _B-band limit at 2 0.79 corresponds to AIb ~ —18.8 at 
2 ~ 0.1. However, the depth of SDSS in B-band is shallower 
than this at 2 > 0.082. Therefore, we are unable to construct 
a volume-limited sample from SDSS matching the depth of 
DEEP2. In principle, we can apply 1/Vmax weighting, but 
it is unnecessary, as the number of K+A galaxies is small, 
which only allows us to make one or two wide bins in lu- 
minosity. A 1/Vmax weighting does not significantly change 
the effective median magnitude of K+As in each bin. In 
addition, the local K+A's environment distribution has no 
detectable luminosity dependence, as shown in fj6l We there- 
fore do not apply any magnitude-dependent weights. 

To maximize the sample from SDSS while making 
it B-band limited, we empirically determined a redshift- 
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Figure 2. H/3 EW vs. /a for the SDSS sample (left) and the DEEP2 sample (right) as defined in Table [T] The solid lines indicate the 
selection criteria adopted for identifying K+A galaxies. The crosses at position [0.8,5] indicate the median 1-cr measurement error for /a 
and H/3 EW in each sample. In the left panel, the K+A galaxies comprise a clearly distinct population and stand out as the horizontal 
spur with nearly zero H/3 EW but relatively large A-star fraction (/a)- The spur is not as obvious in the DEEP2 sample due to the greater 
noise in spectral index measurements and smaller sample. However, as described in i|3.3l we can reliably measure the contamination due 
to scatter in both dimensions and quantify the true K+A fraction. 



dependent restframe B-band luminosity limit. Because 
SDSS is r-band limited, this B-band limit is set by the bluest 
galaxies. The limit we apply is: 

Mb{z) - 51og/i< max(100z^ - 452 - 15.8, -18.8) (3) 

We also select a bright subsample in SDSS with Mb — 
blogh < -19.8, which is the depth of SDSS at 2 = 0.12, the 
upper limit in redshift for our 2 ~ 0.1 sample. This bright 
subsample for SDSS is therefore volume-limited. 

For the environment analysis, we also exclude galax- 
ies within lh~^Mpc (comoving) of a survey edge or hole, 
along with regions on the sky with less than 80% complete- 
ness. With these criteria, we are left with 79716 galaxies, 
amongst them 187 K+A galaxies. The redshift range, ef- 
fective redshift, and magnitude limits adopted for the two 
parent samples are summarized in Table [T] 



3.3 K+A Galaxy Fraction and the Contamination 
Correction 

It is easily seen in Fig.[2]that the K+A population in DEEP2 
is not so obviously a separate population from star-forming 
galaxies as in SDSS. This is largely due to the much larger 
sample size of SDSS, but the greater measurement errors 
in DEEP2 spectral indices also obscure the gap. Therefore, 
to reliably estimate the K+A galaxy fraction, we need to 
do a contamination correction, i.e., estimate the expected 
false-positive rate and the true-negative rate for a given en- 
vironment bin, then subtract or add them, respectively, to 
the raw K+A fraction. 

A robust contamination estimate requires reliable er- 



ror estimates for both spectral indices used. The error on 
each index is first estimated from the inverse variance array 
output from the DEEP2 DEIMOS spectroscopic reduction 
pipeline for each object. However, there can also be system- 
atic errors due to observational conditions, such as incom- 
plete slit coverage, variable seeing conditions, misalignment 
between slit and atmospheric dispersion direction, etc., that 
cause pipeline errors based on measurement noise to be an 
underestimate of the true measurement er ror. To give reli- 
able errors, we follow iBalogh et all l|l999D and use objects 
with repeated observations to scale up the pipeline errors 
such that repeated measurements on the same object give 
consistent results. The errors on H/3 EW and /a are scaled 
up by factors of 1.81 and 1.61, respectively. The details of 
this rescaling can be found in Yan et al. in prep. Similarly, 
the measurement errors on galaxies in SDSS are also scaled 
according to the estimates from repeated observations. 

With reliable error estimates for H/3 EW and /a, we can 
predict how much contamination these errors (presumed to 
be Gaussian) will cause, assuming the true distribution is the 
same as the observed distribution. This is done by Monte- 
Carlo simulations. We add Gaussian noise to the observed 
H/3 EW and /a values, according to the error estimates of 
each point. We then calculate, for each point, the likelihood 
it falls into the K+A definition box. This likelihood is also 
the probability for each galaxy being a true K+A galaxy. 
The difference between the sum of the likelihoods and the 
raw fraction of K+A galaxies gives the contamination cor- 
rection. The fact that the observed distribution differs from 
the true distribution (due to measurement errors) causes 
this contamination estimate to be imperfect. However, for 
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Table 1. Sample definition summary 



Name 


Redshift 


Median Z 


Mag Range 


Median 


Sample 


K+A 


K+A fraction 


Contamination 




Range 




(Af = Afs - 5 log h) 


Mb -5 log h 


Size 


Sample 




Correction 


SDSS-whole 


0.07-0.12 


0.09 


M < -18.8 


-19.73 


79716 


187 


0.23 ± 0.02% 


0.05% 


SDSS-bright 


0.07-0.12 


0.09 


M < -19.8 


-20.16 


35717 


93 


0.26 ± 0.03% 


0.02% 


SDSS-faint 


0.07-0.12 


0.09 


-19.8 < M < -18.8 


-19.42 


43999 


94 


0.21 ± 0.02% 


0.07% 


DEEP2-whole 


0.72-0.88 


0.79 


M < -19.7 


-20.36 


2649 


34 


1.3 + 0.2% 


0.9% 


DEEP2-bright 


0.72-0.88 


0.79 


M < -20.7 


-21.02 


783 


17 


2.3 + 0.5% 


1.1% 


DEEP2-faint 


0.72-0.88 


0.79 


-20.7 < M < -19.7 


-20.12 


1866 


17 


0.9 + 0.25% 


0.8% 



distributions with small second derivatives (as here), this 
effect is modest and may be safely ignored. 



3.4 Environment Measures 

The local galaxy environment has been estimated for each 
galaxy in the DEEP2 and SDSS samples that i s at least lh~ ^ 
Mpc (comoving) away from a survey edge (see 
20061. 12OO7I for measurements in DEEP2 and 
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20081 in SDSS). The environment indicator employed is the 
projected 3'^'^-nearest neighbor surface density (E3), which 
is derived from the projected distance to the 3'^'^-nearest- 
neighbor, I?p,3 , through E3 = 3/(7rDp 3). The measurement 
is made for objects within a velocity window of ±1000km/s 
along the line-of-sight in order to exclude foreground and 
background galaxies. For each galaxy, the surface density is 
divided by the median E3 of galaxies at that redshift to cor- 
rect for the redshift dependence of the sampling rate in both 
the DEEP2 and SDSS surveys. This converts the E3 values 
into measures of overdensity relative to the median density 
(given by the notation 1 + ^3 here). D etails of the e n viron- 
ment measurements can be found in ICooper et all l|2006l . 
I2OO7I . l2008l ). The environment indicator employed here for 
the SDSS sample is differ ent from th e indicators used by 
iBlanton et all l|2005l ) and iHogg et al ] (|2003l . I2OO6I ). These 
differences are not very significant as the rank order of en- 
vironment among galaxies should be similar, a nd environ- 
ments axe measured on similar physical scales in lHogg et al.l 
l|2006l ) and here. For detaile d comparisons of a f ew different 
environment indicators, see ICooper et all (|2005l ). 

We note that this environment indicator measures dif- 
ferent scales in different environments, which also vary with 
redshift. At 2; ~ 0.8, simulations show that, at log(l + 53) > 
0.5, this indicator can be roughly translated as the host- 
ing halo mass; at log(l + ^3) < 0.5, this indicator measures 
the inter-halo environment Typically it measures local den- 
sity on a 2Mpc scale. Although SDSS has a higher sampling 
density than DEEP2, our environment indicator measures 
a slightly larger scale in SDSS than in DEEP2: the median 
Dp,3 is ~ 1.3h~^ Mpc (comoving) in the SDSS sample and 
~ 1.05h~^ Mpc (comoving) in the DEEP2 sample. It is im- 
portant to remember that this environment indicator is nor- 
malized relative to the median overdensity of all galaxies at 
each redshift, therefore is a relative measure of environment. 
This is important when we consider environment evolution 
with redshift. 



4 COLOR-MAGNITUDE DISTRIBUTION OF 
POST-QUENCHING GALAXIES 

Figure |3] shows the colour-magnitude distributions of all 
K+A galaxies in SDSS and in DEEP2 selected by the meth- 
ods described above in i]3.1l (circles) , while the distribution 
of the parent sample is plotted in grey countours for SDSS 

and in grey points for DEEP2. 

As seen in numerous studies (e.g. iBell et al.l |2004| : 
IWillmer et "aLll2006l ). the restframe colour distribution of all 
galaxies is bi-modal at both 2: ~ 0.1 and 2 '-^ 0.8, with a tight 
red galaxy sequence and a diffuse 'cloud' of blue galaxies. For 
DE EP2, we use the same magnitude-dependent colour cut 
as in lWillmer et al.l (|2006l ) to separate red and blue galaxies, 
as indicated by the tilted solid line. 



U -B = -0.032(Mb - 5 log h) + 0.343 



(4) 



One might notice that the colour-magnitude distribution of 
SDSS galaxies is shifted redward relative to that of DEEP2. 
This is due to the aging of the stellar population and the 
decrease in the star formation rate from z ~ 0.8 to 2 ~ 0.1 
iBlantonlbOO^ I. The red sequence is also much more densely 
populated at lower 2. This echoes the observed g rowth of this 
popu l ation from lumino sity function studies (jPaber et al.l 
I2007I : Ib rown et al.ll2007l ). Since galaxies are redder at low 
2, we adopt a redder colour cut for SDSS, which is 0.14 mag 
redder than that for DEEP2, the same as that adopted by 
ICooper et al.l (|2008l ). This cut separates the red and blue 
galaxies at low z reasonably well. 

The K+A galaxies are highlighted as circles on these 
plots. The size of each circle indicates the corresponding 
galaxy's probability of being a true K+A (see ij3.3|l . The 
magnitude limit for each sample is indicated by the vertical 
dashed line in each panel. 

Despite the shifts in colour of the red sequence and the 
blue cloud from DEEP2 to SDSS, the median U - B colour 
of the K+A galaxies remains largely unchanged. This is ex- 
pected as they are selected in the same manner, through 
decomposition of the stellar spectra, thus having the same 
light- weighted mean stellar age. This further demonstrates 
the equivalence between the two K+A samples selected from 
the two surveys. In DEEP2, K+A galaxies mostly occupy 
the red side of the "green valley" and the blue side of the 
red sequence. In SDSS, they are mainly in the "green valley" , 
bluer than most of the red sequence galaxies. This difference 
in their positions relative to the red sequence is entirely due 
to the aging of stellar populations on the red sequence from 
2 ~ 0.8 to 2 ~ 0.1. 
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Figure 3. This plot shows the colour-magnitude distributions of SDSS (loft panel) and DEEP2 (right panel) K+A galaxies. The 
background gray contours or points represent all galaxies within the same redshift ranges with a successful redshift measurement. The 
gray points represent all galaxies within the same redshift range in DEEP2 with a successful redshift measurement. K+A galaxies are 
highlighted as circles with sizes proportional to the degree of confidence of their identification. In each panel, the vertical dashed line 
indicates the absolute magnitude limit adopted in each sample in all following analyses; the tilted solid line indicates the red/blue 
demarcation in each sample, which is 0.14 magnitude redder at the SDSS redshifts (z ~ 0.1) than at the DEEP2 redshifts (z ~ 0.8). For 
the DEEP2 sample at Mb — 51ogh fainter than -20.4, we apply an additional cut, S/N per pixel > 1, before making the K-l-A selection; 
these fainter galaxies are plotted here, but not included in our K+A sample in the analyses below. 



5 ENVIRONMENTS OF POST-QUENCHING 
GALAXIES 

In this section, we investigate the environment distribution 
of K+A galaxies. We face several difficulties. First, the envi- 
ronment measurements have large errors, predominantly due 
to the sparse sampling of the matter distribution provided 
by a galaxy survey. Environment differences between sam- 
ples can only be detected when the sample sizes are large. 
Secondly, our K+A sample size, although being the largest 
available at this redshift, is small for statistical purposes. 
Thirdly, all statistical presentations and tests have limita- 
tions, especially given noisy measurements and small sample 
sizes. To interpret the results, we have to keep these facts 
in mind and combine multiple methods to reach a conclu- 
sion. Therefore, we investigate the environment distribution 
of K+A galaxies in several ways. 

5.1 SDSS 

Before we look at the higher-z sample from DEEP2, we first 
investigate the environments of K+As in the SDSS sample, 
which gives a larger K+A sample and better statistics. 

Our first approach is to compare the K+A fractions 
among galaxy subsamples binned by their local overdensity. 
This is similar to the approach in most previous studies, 
which compared the fractions of K+As between two classes 
of environment: clusters vs. field. 

In Fig. |4l we present the K+A fractions relative to all 
galajcies as a function of environment in the SDSS sample. 



The environment bins are chosen so that each one contains 
an equal number of galaxies (1/5 of all galaxies in this case). 

Since the measurement of a fraction has a binomial dis- 
tribution, the error bars are computed according to the vari- 
ance of a binomial distribution. For small fractions, as for 
K+A galaxies, this is very close to Poisson error. These er- 
rors are underestimated because all the bins are correlated 
with each other due to environment uncertainties. The er- 
rors on the contamination correction for K+A samples is 
assumed to be small, since it is based on empirical error 
estimates on spectral measurements with repeated observa- 
tions. 

As shown in Fig. |31 the blue galaxy fraction decreases 
with increasing local density, while red galaxy fraction in- 
creases with local density, as expected from the color-density 
relation. The K+A fraction shows a similar downward trend 
as blue galaxies. 

In Fig. O we present the environment distributions for 
red, blue, and K+A galaxies using histograms and cumula- 
tive distributions. These histograms first show the level of 
spread in local density of these subsamples, which includes 
both intrinsic spread and measurement error. As expected, 
the blue galaxies have a distribution that is centered at a 
lower overdensity than red galaxies, and the K+A galaxies 
have a distribution that is similar. The same thing can be 
seen in the cumulative distributions: blue galaxies occupy 
lower densities than red galaxies, and the K+As' cumula- 
tive distribution is similar to that of blue galaxies, with only 
small differences at the tails of the distribution (of low sta- 
tistical significance). 
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Figure 4. The fractions of K+A galaxies, blue galaxies, and red 
galaxies as a function of environment for the SDSS sample. The 
raw K+A fraction is shown as the dotted line with error bars. 
The K+A fraction after applying the contamination corrections 
described in i|3.3l is shown as the solid line. The horizontal bars 
indicate the 90-percentile ranges of environment for each bin. The 
blue dashed line and the red dash-dotted line, which are plotted 
with respect to the vertical axis on the right, show the blue galaxy 
fraction and the red galaxy fraction, respectively, in the sample 
as a function of environment. The separation between blue and 
red galaxies i s defin ed here using the cut (in U — B) applied by 
ICooper et al] l |2008t ). which is 0.14 mag redder than the DEEP2 
cut to account for the colour evolution of galaxies. K+A galaxies 
in SDSS show a similar downward trend as blue galaxies. However, 
as discussed later, the K+A fraction trend in this plot depends 
sensitively on the sample used for the denominator, making it 
difficult to compare with other samples. The significance is also 
overestimated due correlations between the bins. 



It is quite obvious that the difFerence between the en- 
vironment distributions of blue galaxies and red galaxies is 
very significant, as is the difference between K+A's and red 
galaxies. However, is the small difference between K+As and 
blue galaxies significant? To address this requires a statisti- 
cal test. A commonly-used non-parametric (i.e., independent 
of Gaussian assumption) test is the Kolmogorov-Smirnov 
test (hereafter K-S test), which measures the maximum dis- 
tance between the two cumulative distributions. Under the 
null hypothesis that the two samples tested are drawn from 
the same population, a larger maximum distance is less likely 
to occur. Throughout this paper, we adopt a threshold sig- 
nificance level of 0.05: i.e., we rule out the null hypothesis 
only when the probability of obtaining a statistic as large 
as that measured or larger under the null hypothesis, the 
p-value, is less than 0.05. In Table [2J we show the K-S test 
p-value results for each subsample. 

One drawback of the K-S test is that it is rather 
insensitive to the distribution tails. A powerful alter- 
native, the Anderson-Da rl ing test (he r eafter A-D test, 
lAnderson fc Darling! Il954l : iPettittI 1 19761 : [ Sinclair fc Spurij 



Figure 5. Shown in the upper panel are the local density distri- 
butions of K+A galaxies (black, solid line), blue galaxies (blue, 
dotted), and red galaxies (red, dashed) in the SDSS sample. The 
histograms for the K+A galaxies are magnified for ease of com- 
parison. The lower panel shows the corresponding cumulative dis- 
tribution for each population. In SDSS, the K+A galaxies show 
a very similar environment distribution as blue galaxies, and is 
significantly different from red galaxies. 



'1988^ is more sensitive to differences on the tails of the dis- 
tribution. In this case, the test statistic is a weighted inte- 
gration of the squared differences between the two cumu- 
lative distributions, with more weight placed on differences 
in the tails of the distributions. In addition, we also apply 
the Mann- Whitney U test (also called the Mann-Whitney- 
Wilcoxon, or Wilcoxon r ank-sum test, hereafter MWW test, 
iMann fc WhitnevI [l947l ) . which is more sensitive to differ- 
ences than the K-S test when sample sizes are small. It uses 
the relative ranks of measurements in a combined, sorted list 
to test the hypothesis that two samples might come from the 
same population. We present the results for all these tests in 
Table[21 Under all of these tests, the null hypothesis that red 
galaxies and K+A galaxies are drawn from the same distri- 
bution is rejected with a significance level of 10~* or better. 
The null hypothesis that the blue galaxies and K+A galaxies 
are drawn from the same distribution cannot be rejected at 
the 0.05 significance level, however. Therefore, we conclude 
that K+A galax;ies at « ~ 0.1 have an environment distribu- 
tion that is indistinguishable from blue galaxies, given our 
samples. 



5.2 Comparison of Methodologies 

All the methods presented above - fraction vs. environment 
plots, density distributions, cumulative distributions, and 
statistical tests - lead to the same conclusion. The K+A 
galaxies at z ~ 0.1 in SDSS have an environment distribu- 
tion which is distinct from that of red galaxies, but similar 
to that of blue galaxies. 
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Table 2. Statistical Tests of Differences Between Samples. This table lists the statistical tests results on the differences between the two 
samples in each pair. The values listed are the probabilities of obtaining a test statistic as large as that measured or larger under the null 
hypothesis that the two samples are drawn from the same population. This probability is called the p- value for most statistical tests. We 
use a threshold signficance level of 0.05, i.e., we rule out the null hypothesis only when the p- value is less than 0.05. 



Sample 


subsamplcs 


Kolmogorov-Smirnov 


Anderson- Darling 


Mann- Whitney U test 


SDSS 


red vs. blue 


< 10-5 


< 10-5 


< 10-5 




red vs. K+A 


< 10-5 


< 10-5 


< 10-5 




blue vs. K+A 


0.419 


0.306 


0.246 


DEEP2 


red vs. blue 


< 10-5 


< 10-5 


< 10-5 




red vs. K+A 


0.739 


0.912 


0.937 




(weighted) 


0.980 


0.946 






blue vs. K+A 


0.166 


0.120 


0.080 




(weighted) 


0.116 


0.049 





The first method, the fraction vs. environment plot, is 
the easiest to read but the hardest to interpret and compare 
between different samples. First, it is very easy to overesti- 
mate the significance of the trend presented in such a plot, 
because the fractions in all the bins are correlated. For ex- 
ample, in Fig. |4j it looks quite convincing that the K+A 
fraction drops consistently, except for the second to last bin. 
In fact, the middle three bins are basically consistent with 
being fiat. Their relative heights change significantly if we 
shift the bins slightly. The most significant variation in the 
K+A fractions is the difference between the leftmost bin and 
the rightmost bin: the latter is significantly lower. 

Secondly, the slope of the trend sensitively depends on 
the relative mix of blue and red galaxies in the sample we 
are taking a fraction of. If we increase the overall fraction 
of red galaxies without changing the K+A population - for 
instance, by switching to an 7-band selected sample - we 
will include more galaxies in overdense environments. Un- 
der the same equal-number binning, the bin boundaries will 
shift to higher densities. Even if the K+A population in- 
cluded does not change, the slope of the K+A fraction will 
still become steeper, i.e., decrease faster towards denser en- 
vironments. The lesson is that the slope of the fraction vs. 
environment trend is sensitive to the parent sample, which 
means we cannot blindly compare the K+A fraction slopes 
between different samples, such as ones defined using dif- 
ferent limiting band, using different luminosity cuts, or at 
different redshifts. We will see an example of this when we 
investigate luminosity-dependence below. 

The second presentation method used is the density dis- 
tribution histogram. The histograms provide much more de- 
tail than the fraction plot, but for small samples like K+A 
galaxies, the number of bins has to be relatively small, thus 
limiting the power of the data. Cumulative distributions give 
a much more sensitive presentation: it does not require bin- 
ning, and thus can take full advantage of the data while mak- 
ing the presentation as continuous as possible. We therefore 
consider this presentation the most accurate and informa- 
tive. 

In fact, many of the non-parametric statistical tests are 
based on comparisons between empirical cumulative distri- 
bution functions (CDF). They provide quantitative tests of 
the significance of the differences between different samples. 
They also have limitations, however. Different tests compare 
different aspects of the data, thus they do not necessarily 



give the same answer. The test giving the smallest p- value 
is not necessarily the most powerful. This will be clear in 
the following analysis for DEEP2 K+As. 

Comparing all these methods, we consider the cumula- 
tive distribution the most accurate and informative way of 
presentation, and we base our judgements on the statistical 
tests for our conclusions. 

Now we move on to investigate the environments of 
K+A galaxies in the DEEP2 sample, which is a smaller sam- 
ple and hence more noisy. 

5.3 DEEP2 

As in >i5.1l we present the K+A environment distribution in 
DEEP2 to blue and red galaxies with three different meth- 
ods. Figure [S] shows the K+A fraction as a function of envi- 
ronment, along with that for blue galaxies and red galaxies. 
The blue fraction decreases toward denser environments and 
the red fraction increases, as expected. The K+A fraction 
shows a slight preference for the central environment bin (al- 
though the contamination correction is large for the DEEP2 
sample, the trend does not change with or without the cor- 
rection). This indicates its environment distribution is in 
between that of the blue and red samples. However, due to 
the correlation among the three bins and the small sample 
involved, the K+A trend presented by this plot is not very 
robust. We turn to the density distributions and cumulative 
distributions for deeper analysis. 

Figure[7]shows the density distributions and cumulative 
distributions for red, blue, and K+A galaxies in the sample. 
For K+A galaxies, we make two sets of density distribution 
and cumulative distribution: one set is based on pure number 
counts, the other set is weighted by the probability of each 
K+A being a true K+A galaxy. We do this because the 
spectral measurements in DEEP2 have larger errors than 
those in SDSS and the sample size is small. Thus, giving 
more secure K+A galaxies higher weight can make the result 
more robust. 

As seen in this figure, K+A galaxies show a density dis- 
tribution skewed towards overdense regions, more similar to 
red galaxies than blue. The cumulative distributions basi- 
cally tracks that for red galaxies, especially when weighted 
by K+A probability. The statistical tests also show that the 
environment distribution of K+A galaxies is more similar to 
red galaxies than to blue galaxies, as seen in Table [2] K+A 
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Figure 6. The fractions of K+A galaxies, blue galaxies and red 
galaxies as a function of environment for the DEEP2 sample. All 
symbols follow the same conventions as in Fig. H| The separa- 
tion b etween blue and red galaxies is defined as in lwillmer et aO 
||2006| ). The error bars for the blue and the red fractions are the 
same since they sum to 1. K+A galaxies appear to have a slight 
preference for the central environment bin. However, the signifi- 
cance of this trend is plagued by small number statistics and is 
less robust. 



Figure 7. Shown in the upper panel are the local density distri- 
butions of K+A galaxies (black, solid line), blue galaxies (blue, 
dotted), and red galaxies (red, dashed) in the DEEP2 sample. 
The histogram for the K+A galaxies is magnified for ease of com- 
parison. The lower panel shows the corresponding cumulative dis- 
tribution for each population. In DEEP2, the K+A galaxies' envi- 
ronment distribution is in between those of red and blue galaxies, 
perhaps more similar to blue galaxies than to red galaxies. 



galaxies are consistent with being drawn from tlie same pop- 
ulation as red galaxies in all tests, regardless of weighting. In 
the case for blue galaxies, the probabilities are approaching 
the threshold for rejecting the null hypothesis. However, the 
differences do not reach our significance threshold (p < 0.05) 
except for the weighted A-D test (p=0.049). Thus, we cannot 
safely conclude that K+A galaxies at z ^ 0.8 have differ- 
ent environment distribution as blue galaxies. We can only 
conclude that they are possibly more similar to red galaxies 
than to blue galaxies. The difference become more signifi- 
cant when we weight more secure K+As more. 

Small sample size and large intrinsic error in the envi- 
ronment measures both affect the significance of the results. 
To improve the significance, larger samples and/or more ac- 
curate environment measures are required. 

5.4 K+As in Groups vs. Fields in DEEP2 

Another way to characterize the local environment of galax- 
ies is to classify galaxies in one of two classes: cluster/group 
members and isolated field galaxies. He r e we u se the DEEP2 
group catalog produced bv lGerke et al.l (|2005l ) to investigate 
the K+A fractions in the DEEP2 groups and in the field. 
We note here that few of our galaxies are in large clusters, 
because the volume of DEEP2 is not large enough and there 
are fewer clusters at 2 ~ 1 than z ~ 0. 

The group catalog used here i s an e xpanded version 
of the one presented in lGerke et al.l (|2005l ). using the com- 
plete DEEP2 datasets rather than the first 25% data. The 



groups are identified using the Voronoi-Delaunay Method 
(|Marinoni et al.ll2002l ). For details of the method, along with 
assessments of the co mpleteness and pu rity of the result- 
ing group catalog, see lGerke et al.l l)2005l l. We calculate the 
K+A fraction in galaxy groups with estimated velocity dis- 
persion of (T > 200km/s, which corresponds to a minimum 
dark matter halo mass of ~ 6 x 10^^ Mq (|Coil et al.l 
I2OO6I ). Out of 671 galaxies with Mb — 51og/i brighter than 
— 19.7 in 329 groups within the redshift range considered, the 
raw K+A fraction is 1.0 ± 0.4%. Among 1613 field galaxies 
outside groups, the raw K+A fraction is 1.5±0.3%. Applying 
contamination corrections, the true K+A fractions in groups 
and in the field are 0.1 ± 0.4% and 0.7 ± 0.3%, respectively. 
We also list the K+A fractions in groups by their richness in 
Table O We reach same conclusions regardless of the varia- 
tions in the group sample: the fraction of K+As among field 
galaxies is slightly higher than that among group galaxies, 
at 1-a level. The significance of the result is limited by the 
small number of K+As. The contamination corrections are 
likely overestimates, which makes the net fraction become 
exceedingly small or even negative. 

5.5 Comparison with Previous Studies 

Numerous authors have studied the environment of K+A 
galaxies in the past. Here we review them and compare their 
results with ours. As we mentioned before, the environment 
preference of K+A galaxixes depends on the methodolgy 
employed, the parent sample used, and the selection critieria 
of K+As. In addition, different studies often adopted differ- 



12 R. Yan et al. 

Table 3. K+A Fractions in Groups vs. Fields in DEEP2 



Description 




Raw K+A Fraction 


Corrected K+A Fraction 


No. of Galaxies 


No. of Groups 


Groups w/ Vdis 


p > 200km/s 


1.04 ± 0.39% 


0.13% 


671 


329 


Groups w/ 4 or 


more members 


0.95 ± 0.36% 


-0.13% 


737 


313 


Groups w/ 2 or 


3 members 


0.87 ± 0.35% 


-0.20% 


690 


486 


Field galaxies 




1.49 ± 0.30% 


0.70% 


1613 





ent environment indicators. All these make it difficult to do 
a fair comparison. Nonetheless, we review the results from 
a few previous studies here. 



5.5.1 Low Redshift 

At low redshift, nearly all studies reached the same conclu- 
sion as ours: K+As at low-redshift have a similar environ- 
ment distribution as blue galaxies, which are predominantly 
fou nd in the field. 

IZabludoff et al.l (|l996l ) first did this study at low red- 
shift in the Las Campanas Redshift Survey. With a K+A se- 
lection based on strong Balmer absorption lines and lack of 
[O II] emission, they found that at least 75% of K+A galax- 
ies are found in the field, not in or near clusters. iBlake et al.l 
l|2004l ) used the 2dF Galaxy Redshift Survey and adopted a 
similar selection criteria to study the environment of K+A 
galaxies. They found the low-redshift K+A galaxies are lo- 
cated predominantly in the field. These two studies both 
selected K+A based on the lack of [O II] emission, which 
would make the sample incomplete for K+As hosting a low- 

luminos ity AGN. 

Quintero et al.l l|2004 ) studied the K+A environment in 
SDSS and used a selection method very similar to ours|f| 
They found that the mean overdensity of K+As on the 1 
and 8 Mpc/h scales are similar to the mean of spiral galaxies, 
most of which live in the field, and is lower than the mean of 
bulge-dominated galaxies, many of which occup y high den- 
sity en vironment. With the same K+A sample, Hogg et al.l 
l|2006l ) employed three different environment indicators: the 
number density in 8 Mpc/h comoving spheres, the transverse 
distance to the nearest Virgo-like cluster, and the transverse 
distance to the nearest luminous neighbour. They found that 
K+A galaxies at low z have a very similar environment dis- 
tribution to star-forming galaxies, with only small, statis- 
tically insignificant differences. We confirm these results at 
low z. 

iGotol (|2005l ') selected K+A galaxies in SDSS with strong 
H5 absorption but without either Hq or [O II]. The addi- 
tional requirement of not having [O II] will make the sam- 
ple more restrictive than what we define here. With this 
seleciton, they found the environment of K+As at 0.5, 1.5 
and 8Mpc scales are consistent with that of field galaxies, 
and inconsistent with that of cluster g alaxies. With basically 
the same sample. ISalogh et al.l l|2005h used the E5, derived 
from the projected distance to the fifth-nearest neighbour 
brighter than Mr = —20, to characterize the environment. 
They concluded that K+A galaxies reside in environment 



^ In place of our cut in H/3 EW and f ^ space. [Hogg et al"] 1I2OO6I ) 
used Ha EW and A/K ratio to select K+As. 



typical of normal SDSS galaxies and that are inconsistent 
with overdense regions like clusters. 

5.5.2 High Redshift (0.3 < z < 0.9) 

At redshift greater than 0.3, the environment of K+A galax- 
ies is an unsettled question. But all the high-z studies done 
in the past used [O II] as the star formation indicator to 
select K+A gala xies. As mentioned in H3.1l and discussed in 
lYan et al] (|2006l ). [O II]-based selection methods suffer from 
incompleteness caused by AGN line emission and contami- 
nation from dusty star-forming galaxies. This can be partic- 
ularly bad in clusters since [O II] can be str ong in LINERs 
most of which r eside in early- type galaxies ||Ho et al.lll997l : 
I Yan et all [200^ ). Besides this selection difference, to com- 
pare our results with previous work also requires matching 
the all-galaxy sample. This is critical and was often ignored. 
All these difficulties prevent us from doing a fair quantita- 
tive comparison. Nonetheless, we list the results from a few 
repr esentative work s here . 

iDressler etlii] l|l999l ). iPoggianti et al.l l|l999l ) studied 
the K+As in 10 rich galaxy clusters in the range 0.37 < 
z < 0.56. They found that the K+A fraction in clusters 

18%) is s ignific antly higher than that in the field (~ 2%). 
iTran et al.1 (|2003h studied K+As in 3 clusters at z = 0.33, 
0.58, and 0.83 and deter mined that t he K+ A fractions in 
them are between 7-13%. iTran et al.l (|2004l l quantified the 
K+A fraction in the field sample from the same survey to 
be 2.7 ± 1.1%, thus they concluded that K+A fraction is 
significantly h igher in clusters than in the field. However, 
iBalogh et all |l99i) used CNOCl survey and found simi- 
lar K+A fractions between cluster and field environments. 
The difference from these results are likely due to differ- 
ences in cluster selection (optical selection vs. X-ray selec- 
tion) and in galaxy sample selection (morphological selec- 
tion vs. magnitude limited s ample), which are discussed in 
detail in lBalogh eTsU l| 19991 ). Although we do not have any 
Coma- like rich cluster covered in our sample, we do find the 
fraction of K+As in groups is similar to, perhaps slightly 
smaller than, the fraction o utside groups . This is roughly 
consistent with the result of iBalogh et al.l l|l999l ). However, 
we cannot make finer comparisons due to selection differ- 
ences and lac k of large cluster s in ou r sample. 

Recently, IPoggianti et al.l (|2009l ) studied K+A galaxy 
fractions in clusters, groups, a nd fields using data from the 
ESQ Distant Cluster Survey l| White et al.l l2005l ). Using a 
K+A selection based on [O II] emission and H5 absorption, 
they found that the K+A fraction is higher in clusters and a 
subset of groups with low fraction of [O II] emitters, than in 
the field. Albeit the selection differences, the K+A fraction 
trend they found are actually consistent with our findings. 
If we count both their high-[0 II] and low-[0 II] groups, the 
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K+A fraction among all group galaxies (5 ± 3%) is com- 
parable to the fracti on in the field (6 ± 3%) (Table 4 of 
IPoggianti eral]l2009l ). This is consistent with the compara- 
ble raw K+A fractions we find among DEEP2 group galaxies 
(1.0 ±0 4%) and field galaxies (1.5 ±0.3%). The larger frac- 
tions in lPoggianti et al] l|2009l ) than ours are probably due to 
their more inclusive selection criteria and different all-galaxy 
sample definition. We are unable to compare to their K+A 
fractions in clusters due to the paucity of rich clusters in the 
DEEP2 volume. The Red-Sequence Cluster Survey (RCS), 
with its spectroscopic follow-up prog ram, might prov ide an 
ideal parent sample for such studies (|Yee et al.ll2007^ . 



6 TESTING FOR LUMINOSITY 

DEPENDENCE OF THE ENVIRONMENT 
DISTRIBUTION 

In this section, we test for luminosity dependence in the 
K+A galaxies' environment distribution. We also illustrate 
some serious drawbacks of relying on measurements of K+A 
fraction as a function of environment. 

In Fig. |8l we show the K+A fractions relative to all 
galaxies as a function of environment for two SDSS sam- 
ples binned by luminosity: the SDSS-bright sample {Mb < 
-19.8) and the SDSS-faint sample (-19.8 < Mb < -18.8). 
Comparing the two panels, the blue fraction trend and the 
red fraction trend do not change greatly with luminosity, 
except for the overall fractions of each type. However, the 
K+A fraction vs. local density trends appear quite diff'erent 
in the two samples. For the bright sample, K+A fraction 
decreases with increasing local density, mimicking the envi- 
ronmental trend of blue galaxies. For the faint sample, K+A 
fraction is fairly flat across all environment bins, unlike ei- 
ther the red galaxies or the blue galaxies. Does this mean 
there is a luminosity dependence in the K+A environment 
distribution? The answer is no, for interesting reasons. 

Contrary to one's intuition, the apparent luminosity de- 
pendence of the K+A fraction trend does not necessarily 
reflect a change in K+A galaxies' environment distribution 
with luminosity. In fact, it is largely due to the change in the 
environment distribution of the parent sample used to calcu- 
late fractions. Fraction measurements are dependent on the 
environment distribution of the all-galaxy sample used, not 
just the K+As. In the faint sample, there is a higher fraction 
of blue galaxies and a lower fraction of red galaxies than in 
the bright sample. The median overdensity of the sample 
decreases slightly; since we require that each bin has a equal 
number of galaxies, the bin limits shift slightly to the left. 
This slight shift can change the apparent trend signiflcantly. 

We demonstrate this by replotting the panel for the 
faint sample, but this time replace the K+A galaxies in the 
faint sample by those in the bright sample. Therefore, com- 
paring to the bright-sample plot (left panel of Fig. [8l , the 
K+A galaxies are the same, only the parent sample changed. 
This is shown in Fig.|9] Although the K+A galaxies included 
are the same, the K+A fraction trend displayed is different, 
appearing flat except for the last bin. This demonstrates 
that the trend displayed in a fraction vs. environment plot 
is sensitive to the comparison sample. Comparing this with 
the original SDSS-faint sample (right panel of Fig. [Sj, the 
only difference that may be significant is the K+A fraction 
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Figure 9. This plot demonstrates the sensitivity of the K+A frac- 
tion slope on the parent sample. It is made with the SDSS-faint 
sample, but with the K+A galaxies in it replaced with those in 
the SDSS-bright sample. Comparing this plot with the left panel 
of Fig. [S] which have the same K+A sample but different parent 
sample, we see the K+A fraction slope changed significantly. 



in the last bin. However, we cannot be sure since the er- 
ror bars are underestimates due to correlations amongst the 
bins. 

In Fig. 1101 we show the cumulative distributions of 
environments of the K+As, blue galaxies, and red galax- 
ies in the SDSS-bright and SDSS-faint samples. The two 
K+A CDFs track each other closely in the low-density 
regime. There is an excess in the bright K+A sample at 
log(l + 53) = —0.2; above that, the two curves still follow 
each other in slope, until the faint K+A sample catches up 
around log(l + ^3) = 1.35. The two all-galaxy CDFs show 
very small but statistically significant differences. It is pre- 
cisely this small difference between the two parent samples 
that shifts the binning and makes the slopes of the K+A 
fraction trends appear different. 

The problem, though related to the way in which sam- 
ples are binned, is not solvable by binning in a different 
manner. If, instead of using equal-number binning, we fix 
the binning between the bright sample and the faint sam- 
ple, then the K+A fraction trends can be compared with 
each other. However, the proportion of the parent sample 
in each bin would differ between samples. Thus, any differ- 
ence in the resulting K+A fraction trends may not reflect 
the difference between two K+A samples, but the difference 
between the two parent samples. 

Is there any real luminosity-dependence in the K+A 
environment distribution? Due to the large error in the en- 
vironment measurements, which makes the fraction plots 
useless, we have to rely on statistical tests to answer this 
question. Table. |4] shows the p- values from the same three 
statistical tests used above. The environment distributions 
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Figure 8. The fractions of K+A galaxies, blue galaxies, and red galaxies as a function of environment for the SDSS-bright sample (left) 
and the SDSS-faint sample (right). All symbols follow the same conventions as in Fig.|4] 



Table 4. Statistical Tests of Luminosity Dependence in Environment Distributions 



Sample 


bright vs. faint 


Kolmogorov-Smirnov 


Anderson-Darling 


Mann- Whitney U test 


SDSS 


All 


< 10-5 


< 10-5 


< 10-5 




red 


< 10-5 


< 10-5 


< 10-5 




blue 


4 X 10-5 


< 10-5 


4 X 10-5 




K+A 


0.321 


0.236 


0.333 


DEEP2 


All 


6-5 


< 10-5 


< 10-5 




red 


0.010 


0.0057 


0.0053 




blue 


0.044 


0.0057 


0.0057 




K+A 


0.387 


0.213 


0.823 



of bright K+As and faint K+As are consistent with being 
drawn from the same population. 

We repeat this exercise for the DEEP2 sample. We show 
the fraction plots in Fig. [TT] and the cumulative distribu- 
tion plots in Fig. 1121 Again, the difference in K+A fraction 
trends appear at first to be very significant, but actually is 
not. The two K+A samples are consistent with being drawn 
from the same population based on all statistical tests used 
(see TablelU. However, for the two all-galaxy samples, there 
is significant luminosity dependence in their environment 
distributions. This luminosity dependence of the all-galaxy 
sample is much stronger than that in the SD SS, consistent 
with the results shown in lCooper" et all (|2006| ). 

In summary, with the current samples, we do not find 
a significant luminosity dependence in the environment dis- 
tribution of K+A galaxies at either z ~ 0.1 or 2; ~ 0.8. 

The cautionary lesson we learned here also applies to 
other works looking at K+A fractions as a function of envi- 
ronment. It is vital, when comparing results between differ- 
ent measurements, to ensure that the samples out of which 



K+A fractions are defined are closely matched. Addtionally, 
more care must be taken to assess the significance of these 
results if the environment indicator used has large uncer- 
tainties. These issues may largely explain the differences in 
past K+A environment results in the literature. 



7 DISCUSSION 

The generic quenching picture we assume here is this: a blue, 
star-forming galaxy undergoes an undetermined mechanism 
or mechanisms which may or may not depend on environ- 
ment; the star formation in the galaxy is quenched and the 
galaxy appears as a K+A galaxy; it becomes red in ~ 1 Gyr 
and evolves passively afterwards. 

T.l Redshift Evolution of Environments 

We have seen above that the environment distribution of 
post-quenching galaxies at low-z is similar to blue galax- 
ies and dramatically different from red galaxies, while the 
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Figure 11. The fractions of K+A galaxies, blue galaxies and red galaxies as a function of environment for the DEEP2-bright sample 
(left) and the DEEP2-faint sample (right). 



equivalent population at high-z has an environment distri- 
bution more similar to red galaxies than to blue galaxies. In 
this section, we try to understand this. 

First, we want to remind our readers that the environ- 
ment measure we adopted is a relative environment mea- 
sure. It is computed relative to the median environment of 
all galaxies at each redshift. It is thus different from an ab- 
solute environment measure, such as the dark matter halo 
mass, or the galaxy number density inside a comoving vol- 
ume centered on the galaxy in question. In other words, the 
relative environment is similar to ranking in environment, 
with galaxies in relatively overdense environments having 
higher ranks. This is important for understanding the dif- 
ference in environment distribution between high and low 
redshift. Identical relative environments, i.e., same ranks, at 
different redshifts do not correspond to the same absolute 
environments. 

For the discussion of quenching mechanisms, we shall 
primarily focus on the absolute environment. However, there 
is one unique feature of relative environment. Galaxies are 
less mobile in relative-environment space than in absolute- 
environment space. For example, a galaxy may move into a 
higher density environment (absolute), but its environment 
rank may not change significantly, not more than the un- 
certainty associated with the measurement. This is because, 
statistically, all those galaxies with the same rank have an 
equal probability to have their local densities increased over 
the same time period. Therefore, supposing we can mea- 
sure the environment to arbitrarily high accuracy, then we 
can assume that the environment ranks are fixed statisti- 
cally for all galaxies from z ~ 1 to z ~ 0. In reality, we 
cannot measure environment to arbitrary accuracy, and due 
to the stochasticities involved in this measurement, there 
will be shuffling of ranks among galaxies with time. In ad- 
dition, once the local galaxy density becomes so high that 



we can no longer distinguish differences from redshift-space 
data, for instance, when a group merges into a cluster, the 
ranks hit a ceiling and become meaningless. However, on 
a scale coarser than the measurement uncertainty, the av- 
erage ranks for galaxy populations from different environ- 
ments will be fixed. 

Because environment rank cannot change significantly 
over time, it can be very useful in identifying progenitor- 
descendant relationships. For two populations to have 
a progenitor-descendant relationship, they have to have 
roughly the same environment ranks. 

For the discussion of K-l-A environment evolution, it is 
helpful to visualize the environment evolution of all galaxies 
in a schematic diagram. The left panel of Fig. [13] shows the 
movement of galaxies in ranked environment from z ~ 1 to 
z ~ 0. Overall, environment ranks of galaxies do not change 
significantly over cosmic time. Red galaxies, on average, al- 
ways rank higher in environment than blue galaxies. Our 
results show that K+A galaxies at z ~ 1 are ranked similar 
to red galaxies, but the equivalent population at z ~ has 
much lower ranks and is similar to blue galaxies. This gives 
us two messages. First, K+As at z ~ are a very different 
population from the progenitors of low-z red galaxies. Sec- 
ond, the building-up of the red sequence through the K+A 
phase is happening in relatively overdense (higher rank) en- 
vironments at 2; ~ 1 but in relatively underdense (lower 
rank) environments at z ^ 0. This statement may or may 
not apply to the overall build-up of the red sequence, de- 
pending on whether the K+A phase is the dominant route 
for the production of a red galaxy, and on the environment 
evolution of other routes. Nonetheless, our result at z ~ 1 is 
consistent with the result s from other e nviro nment studies. 
ICooper et al] (120071 ) and iGerke et all l|2007l ) both showed 
that between redshifts 0.7 and 1.3 the red sequence was pref- 
erentially built up in overdense environments. 
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Figure 13. These schematic diagrams show the movements of galaxy environments with redshift. The left panel shows the movements 
in relative environment or ranked environment. The right panel shows the movements in absolute environment (e.g., dark matter halo 
mass, or the local comoving number density of galaxies). Thick red (blue) arrows represent the movements of red (blue) galaxies. K+A 
galaxies at low-z are found at lower rank environments than high-z K+As. See ijTlfor full discussion. Note, the dashed arrows linking 
high-2 K+As with low-z K+As do not represent a progenitor-descendant relationship: they are completely different populations. 



Although high-0 K-|-As have much higher environment 
ranks than low-z K-|-As, they may have the same absolute 
environments. This is because the absolute environments for 
most galaxies have grown with time. In regions with matter 
density above the critical density, the density will grow with 
time, and so will the galaxy number density, if galaxies trace 
the dark matter distribution on large scales. Most galaxies 
are in such regions; thus their absolute environments will 
grow with time. Therefore, the lower environment rank of 
the K-fA population at lower redshift does not necessarily 
indicate a difference in absolute environment from high-z 
K-I-As. This is iUustrated in the right panel of Fig. [I3l The 
absolute environments of K-|-As at z may be the same, 
or less dense than, the absolute environments of K-|-As at 
z 1; but they cannot be much denser than z ~ 1. Other 
than that, we do not have any constraint on this because we 
only have a relative environment measure. 

If we assume K-(-As at all redshifts are quenched in 
the same absolute environments and every red galaxy went 
through this phase, then we can explain the larger envi- 
ronment difference between red galaxies and K-f-As at low z 
than that observed at high z. The average age of stellar pop- 
ulations in red g alaxies is much older on average at z ~ 0. 1 
than at z - 0.8 l|Gallazzi et al.ll2005l : [Schiavon et al.ll2006l ). 
A typical low-z red galaxy would probably have quenched a 
longer amount of time ago, thus its local environment would 
have more time to grow (e.g., the group it belongs to could 
accrete more members and/or merge into a cluster) com- 
pared to a high-z red galaxy. This will lead to a larger dif- 
ference in environment with K-|-As at low-2 than high-z, 
as observed: high-z red galaxies still reside in environments 
similar to those of K-|-As, maybe only slightly denser; while 
low-z red galaxies reside in much denser environments than 
low-z K+As. 

This picture not only explains our results with regard 



to red galaxies, it also has a prediction. Under this scenario, 
the absolute environments of red galaxies will be correlated 
with the time passed since quenching, and hence, approxi- 
mately, the stellar age. If confined to a single redshift epoch, 
the same statement can be applied to relative environments 
as well. This will be investigated in a forthcoming paper 
(Cooper et al. in prep). 



7.2 Quenching Probability Among Blue Galaxies 

The blue galaxy population must include the direct progen- 
itors of future K-l-A galaxies. Comparing the environment 
distribution of K-|-As with that of blue galaxies will reveal 
how quenching may depend on environment and/or inter- 
nal properties of galaxies. This dependence also needs to 
produce the observed evolution in K-l-A environments with 
regard to that of blue galaxies. We discuss these topics in 
this section. 

Only galaxies which possessed a significant (by fraction) 
young stellar population before quenching can appear as a 
K-hA after quenching. If K-fA galaxies are produced from 
a starburst and then quenching associated with a merger 
event, this would be a constraint on the amount of gas avail- 
able for star formation (as a fraction of stellar mass) before 
the merger. Hence, the observed environment distribution of 
K-l-A galaxies will be a result of both any environment de- 
pendence in the quenching mechanism and any environment 
dependence in pre-quenching star formation in galaxies. To 
further complicate matters, the quenching probability could 
depend on galaxy properties which correlate with environ- 
ment, such as stellar mass or SFR, as well as the environment 
itself. 

At z ~ 0.8, we found that the environment distribution 
of K-l-A galaxies is marginally different from that of blue 
galaxies. However, the significance of the difference is lim- 
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Figure 10. Shown in the upper panel are the cumulative density 
distributions of K+A galaxies and all galaxies in the SDSS-bright 
and SDSS-faint samples. Due to the small samples of K+As when 
split in the manner, the difFerenee between the two K+A curves 
are not statistically significant (see Table. |4j: they are consis- 
tent with being drawn from the same population. The difference 
between the two all-galaxy curves, albeit small, is very signifi- 
cant. This is the major cause for the slope difference in the K+A 
fraction trend between the bright and faint sample. The lower 
panel shows the cumulative overdensity distributions for the blue 
galaxies (colored blue, on the left) and red galaxies (colored red, 
on the right) in the SDSS-bright (solid) and SDSS-faint (dashed) 
samples. This nicely illustrate that the environment distribution 
depends on both galaxy color and luminosity, but the dependence 
on the former dominates. 



DEEP2 




Log(l+(53) 

Figure 12. Shown in the upper panel are the cumulative density 
distributions of K+A galaxies and all galaxies in the DEEP2- 
bright and DEEP2-faint samples. Due to the small samples of 
K+As when split by luminosity, the difference between the two 
K+A curves are not statistically significant (see Table. |4j: they 
are consistent with being drawn from the same population. The 
difference between the two all-galaxy curves is highly significant. 
This is a major cause for the difference in the K+A fraction trend 
vs. environment between the bright and faint samples. The lower 
panel shows the cumulative distributions for the blue galaxies 
(colored blue, on the left) and red galaxies (colored red, on the 
right) in the DEEP2-bright (solid) and DEEP2-faint (dashed) 
samples. This nicely illustrates that the environment distribution 
depends on both galaxy color and luminosity, but the dependence 
on the former dominates. 



ited by the small sample size. We must wait for future data 
to confirm this result and improve its significance. In the 
meantime, we will assume that the difference is robust for 
the purpose of discussion. There are three possible explana- 
tions for this difference. 

The difference can be explained as a dependence of 
quenching probability on the local environment. Assuming 
environment is the only factor determining quenching prob- 
ability and all blue galaxies satisfy can appear as a K+A 
after quenching, our results indicate that quenching prob- 
ability at z ~ 0.8 is higher in denser environments, but at 
2 ~ 0.1, it is roughly the same for blue galaxies in any envi- 
ronment (where we here mean absolute, rather than relative, 
overdensity). The trends at the two redshifts can be consis- 
tent only if, by 2 ~ 0.1, nearly all blue galaxies have left the 
low-density regions where quenching is less likely to happen. 
Whether the overall increase in absolute environment for all 
blue galaxies is large enough to make this difference awaits 
simulations to answer. 

The second possibility is that star formation quenching 
does not depend on the local environment, but quenching 
and the K+A prerequisites depend on certain internal prop- 
erties of a galaxy, which selects preferentially blue galaxies 
that are found in overdense environments. For example, if 



quenching requires a sufficiently high stellar mass (e.g., to 
have a massive enough central SMBH to provide the feed- 
back) and forming a K+A requires a high pre-quenching 
SFR, the K+A population will carry the environment foot- 
print of those galaxies satisfying these criteria. At 2 ~ 0.8, 
these galaxies do pr eferentially res i de in overdense environ- 
ments. As shown bv lCooper et al.l (|2006l ). the average over- 
density on the blue cloud at z ~ 1 is a strong function 
of B-band luminosity, which correlates strongly with stel- 
lar m ass, with som e sensitivity to the presence of young 
star s (ICooper et al. ' '20081. In addition. lElbaz et al.l |2007l ) 
and lCooper et al!r ( 2008. ) have both shown that the average 
SFR increases with galaxy overdensity at 2 ~ 1. Therefore, 
K+As at 2 ~ 0.8 will be found in more overdense environ- 
ments than blue galaxies. 

At 2 ~ 0.1, the average environment of galaxies on 
the blue cloud has no or only weak dependence on "'Ri- 
band luminosity (|Hogg et al.ll2003l : lBlanton et al.l2005l ). The 
aver age SFR is als o higher in underdense environments 
( Cooper et al.| [2C)08l). As a consequence, if K+A progenitors 
must be high- SFR blue galaxies, and quenching probabil- 
ity increases modestly with density, we might expect low- 
redshift poststarburst galaxies to have an environment dis- 
tribution similar to blue galaxies. Therefore, the observed 
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environment distribution of K+As relative to blue galaxies 
may indirectly be a reflection of the correlations between 
environment and other galaxy properties. 

This possibility can be distinguished from the first by 
controlling each galaxy property that has an observable cor- 
relation with environment. It is also likely that quenching 
depends on both environment and internal properties of 
galaxies. Future surveys yielding much larger sample sizes 
and higher completeness are required to break the possible 
degeneracies. 

The third possibility is that quenching depends on hav- 
ing a fast and significant increase in local overdensity. Only 
those galaxies that experience a quick change in environ- 
ment (e.g., by falling into a cluster) have their star formation 
quenched. Ram-pressure stripping in clusters is an example 
of a mechanism that has such dependence on environment. 
We would expect K-|-As to be found in rich groups and clus- 
ters in this case. This may occur, but it does not explain the 
field K-|-As. Additionally, it will be difficult to explain the 
low-z result, where the K-fAs have an environment distri- 
bution indistinguishable from blue galaxies. 



8 IMPLICATIONS FOR K+A QUENCHING 
MECHANISMS 

What does the environment distribution of K-fAs tell us 
about their formation? First of all, if a quenching mecha- 
nism has a dependence on environment, the environment we 
should consider is the absolute environment, not the ranked 
environment. Therefore, we have less constraining power on 
this issue as our environment measurement is a relative mea- 
sure. Furthermore, the fact that the ranked environment dis- 
tribution of K-|-As at 2 ~ 0.8 differs from that at 2; ~ 0.1 
does not mean that the quenching mechanisms have to be 
different. In fact, it is possible that K-f As reside in similar 
absolute environments at both redshifts, and are due to the 
same quenching mechanism. 

A variety of quenching mechanisms that can yield 
K-l-A galaxies have been proposed since this population was 
identified. One class of mechanisms operates only in rich 
cluster environments, such as ram-pressu re stripping (e.g. 
iGunn fc GottI Il972l: iBekki fc CouchI |2003| ). g alaxy harass- 
ment ( Moore et al.lll996l ). and strangulation (|Balogh et al.l 
I2OO0I ). Another popular mechanism is through ma j or me rger 
events l|Mihos fc HernauistI [l994l : ISpringel et all |2005|). In 
addit ion, shock heating in massive halos i Birnboim et al.l 
I2OO7I ) may also produce K-l-A galaxies. Below we discuss 
the relevance of each mechanism and our constraints on its 
significance at both z ^ Q and z ^ 1. 



8.1 Cluster-specific mechanisms 

The first mechanism proposed for K-fA formation results 
from interactions of an infalling galaxy with the gravita- 
tional potential or intracluster medium of a rich galaxy clus- 
ter. The hot intracluster medium could shock or compress 
the gas in a galaxy producing a starburst, while also strip- 
ping the gas in its outskirts and other low d ensity parts of 
the ISM, eventually halting star formation (|Gunn fc Got! 
19721 : lNuIsenlll982l: iBothun fc Dressie3ll986l : iBalsara et al 
1994l : lBekki fc CouchI booi i 



Many examples of galaxies undergoing strong ram pres- 
sure stripping have been found in nearby clusters (e.g., 
Virgo, Coma, A1367); asymmetric radio conti nuum mor- 
phology, HI profile, and/or ionized gas t ails (e.g. frwin et al.l 



1987; VoUmer et al. 



McConnachie et al. 



■2000l . |2004 I2OO5I : ICortese et all 12007 : 
2007h have provided evidence for recent 



stripping events. However, the star formation in these galax- 
ies is not always shut down. In some case s, the specific 
SFR is enhanced (CGCG 97-073 and 97-079. iGavazzi et al.l 
Il995h or ev en raised to starbur st levels (235144-260358 in 
Abell 2667. ICortese et al.|[2007l ): some objects show signs 
of moderate SF or have SF rejuvenated from stripped gas 
falhng back (e.g. NGC 4522, which has SF in the inner 
disk and an extraplanar ring, but a K-|-A-like outer disk , as 



found bv lVoUmer et al.|l2000l andlCrowl fc Kennevll2006l : or 
NGC4848. as found bv lVollmer et aLlbOOll ). In two objects, 
however, there is no ongoin g SF and th e systems appear as 
K+A galaxies (NGC 4569 



AbeU 1689, ICortese et al 



Keell 



19961 . and 131124-012040 



20071 ) . All these observations 



suggest that ram-pressure stripping can make K-l-A galaxies 
in some rich groups or clusters, but not in every case. This 
is expected, as the gas in the inne r part of massiv e disks 
is difficult to strip away. Studies bv lPoggianti et ahl (|2004h 
found that K-l-A galaxies in Coma are all low-luminosity ob- 
jects (M„ > —18.5), quite different from the high-luminosity 
K-|-A's studied in this paper {Mb < —19.7). If ram pressure 
stripping is responsible for both low-z and high-z K-|-As, 
it should only be more difficult to strip larger galaxies at 
higher redshifts as cluster masses were lower then. Therefore, 
it is not clear from this evidence that ram-pressure stripping 
could be the dominant K-l-A formation mechanism at high 
z. 

Because of the paucity of rich clusters in the DEEP2 
volume, we are unable to directly address this question at 
high-z. It is possible that cluster-specific mechanisms can 
make K-|-As in clusters, but certainly that is not the domi- 
nant mechanism for K-l-A formation. Regardless, given the 
presence of K-|-As in low-density regions at z ~ 1, and 
the predominance of K-|-As at z ~ 0, cluster-specific mech- 
anisms cannot explain the observed characteristics of the 
K-f A population. 



8.2 Galaxy- Galcixy Mergers 

A second mechanism proposed for K-l-A formation is galaxy- 
galaxy merg ers. Mergers can trigger br ief, intense starbursts 
in galaxies l|Mihos fc HernauistI 1 1994! ) ■ The exhaustion of 
available gas and feedback fro m supernovae winds or AGN 
activity (Springel et all |2005| ) could then halt star forma- 
tion. Many K-|-As have been observed to have morphologi- 
cal features characteristi c of mergers, such as di s turbed mor- 
phologies and tidal tails l|Liu fc Kennicuttll995l ). lYang et al.l 
(2008) recently presented the detailed morphologies of 21 
low-z K-f A galaxies obtained from HST/ACS and WFPC2 
observations. They found that more than half of the K-f-A 
galaxies have disturbed morphology and/or tidal features, 
indicating recent merger events. 

Is the environment distribution observed for K-fA 
galaxies consistent with this mechanism? Merger rate prob- 
ably depends on environment. However, as we mentioned 
earlier, environment may not be the only factor in deter- 
mining the probability a galaxy evolves into a K-fA through 
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merger. For the final product of a merger event to quahfy 
as a K+A, the galaxies involved has to be relatively gas- 
rich. The star formation quenching effectiveness after the 
merger can also depend on the ex istence and mass of a cen- 
tral SMBH (|Springel et al.ll2005l ). which in turn could de- 
pend on the final stellar mass. 

The most favorable environment for galaxy- galaxy 
mergers is in poor groups becau se of the low ve locity disper- 
sion in these systems ( Zabludoff fc MulchaevI [l998 ) . How- 
ever, even poor groups with masses as low as 6 x 10^^ Mq 
in our sample generall y have lo g(l -|- 5^) > 0, reflecting their 
strong clustering |CQiLet_a],' 2 0061 ) : and the K-l-A fraction 
in these poor groups is lower than among galaxies in isola- 
tion. It seems that these facts disfavor the merger scenario. 
Actually, because only gas-rich mergers is capable of pro- 
ducing K-|-As, what we should be looking at is the K-l-A 
rate relative to gas-rich galaxies only, or blue galaxies. 

Relative to blue galaxies, K-|-As at z ~ 0.8 are more 
frequently found in denser environments. This can be con- 
sistent with the merger scenario, although the preference for 
dense environment could partly be due to a quenching re- 
quirement on stellar mass. At 2 ~ 0.1, K-l-A rate relative 
to blue galaxies probably varies little with environments. To 
explain this under the merger scenario, all blue galaxies at 
2 ~ 0.1 have to have equal probabilities to have a major 
merger. This is not impossible, but we need simulations to 
answer this, which is missing in the current literature. 

8.3 Halo Quenching 

A third quenching mechanism that could produce a K-l-A 
galaxy, but has not often been mentioned in K-l-A literature, 
is through shock heating of the halo gas and the inflowing gas 
when a halo grows above a certain threshold in ma s s. Th e 
idea dates back to the early work of lRees fc Ostrikeil (| 19771 ) , 
which concluded that halos with M > 10^^ Mq would 
experience a quasi-static co ntraction phase due to virial 
shocks. A series of studies by iBirnboim fc Dekel (|2003l ) and 
IPekel fc Birnboim ! (2006) using spherical halo gas-accretion 
simulations, and by O cvirk et al,: (^OOS) using N-body and 
hydrodynamic simulations, concluded that when a dark 
matter halo grows above a critical mass (McHt ~ 10^2 Mq), 
shocks will form in these halos, slow down the infalling ma- 
terial and heat it up, thus removi ng the cold gas supply 
for star formation. Recent work by IBirnboim et al.l (|2007l ) 
has noted that the whole process actually involves a se- 
quence of quenching and bursting events due to shock in- 
stabilities at early stages. In cases where shocks form rel- 
atively late - at z < 2.5 - and the subsequent burst is at 
z < 1.4, the final post-burst shock heating is both rapid 
and effective, and it provides long term suppression of cold 
gas accretion onto the central disk of the halo. However, 
the implications of this model for the SFR in the central 
galaxy are not straightforward. The gas accretion rate on 
the central disk provides only an upper limit to the SFR; 
the actual SFR evolution could be more complicated. How- 
ever, the overall burst-foUowed-by-quenching trend should 
remain robust. Therefore, this model provides a mechanism 
for forming K-l-A galaxies. If it is the dominant mechanism, 
K-l-A galaxies will mostly be found in halos that just crossed 
the mass threshold. 

It is possible for this mechanism to explain the environ- 



ment evolution of K-|-As, since the absolute environment for 
K+As may not change signiflcantly from high-z to \ow-z. 
Same halo mass will correspond to different relative envi- 
ronments at different redshifts. This can be tested with a 
simulation that has the correct clustering strength for both 
blue and red galaxies, at both z ~ 0.1 and z ~ 0.8. We defer 
this to future work. 



9 SUMMARY AND CONCLUSIONS 

We have compared the environment distributions of the 
post-quenching/K+A galaxies with those of red and blue 
galaxies at both z ~ 0.1 (from SDSS) and z ~ 0.8 (from 
DEEP2), using uniformly selected samples defined by an 
H/3 EW limit rather than a less-robust [O II] EW limit. We 
have found the following facts. 

(i) At high-z {z ~ 0.8), the environment distribution of 
the K-l-A galaxies is indistinguishable from that of the red 
sequence galaxies, preferring overdense environments, but 
showing modest difference from that of the blue galaxies. 

(ii) At \ow-z {z 0.1), K4-A galaxies are found to in- 
habit environments dramatically different from those of red 
galaxies. Instead, local post- quenching systems are observed 
to have similar environment distribution as blue galaxies, 
preferring underdense environments. K-|-As at low redshift 
have lower rank in the environment distribution of bright 
galaxies than their high redshift counterparts do. 

(iii) We do not find any significant dependence on lu- 
minosity for the K-l-A environment distribution, at both 
z ~ 0.1 and z ~ 0.8. 

(iv) At 2 ~ 0.8, the K-l-A fraction among DEEP2 group 
galaxies is consistent with the fraction among field galaxies 
within 1-CT. 

Based on these facts, we have reached the following con- 
clusions. 

(i) The quenching of star formation and the build-up of 
the red sequence through the K-f A phase is happening in 
relatively overdense environments at 2: ~ 1 but in relatively 
underdense environments at 2 ~ 0. The absolute environ- 
ments (i.e., local mass density) for these processes, however, 
may stay the same with time, as may the quenching mech- 
anisms. This can be true while the relative quenching envi- 
ronment decreases because the mean absolute overdensities 
of all galaxies has to increase with time due to gravity. 

(ii) At 2; ~ 0.8, the typical environment for quenching has 
not changed significantly from the time L* red galaxies at 
2 ~ 0.8 were quenched, suggesting K-|-As at 2 ~ 0.8 may be 
similar to the progenitors of red galaxies found at the same 
redshift. In contrast, the low-2 K-|-As, having much lower 
rank in environment distribution than red galaxies, have to 
be a totally different population from the progenitors of low- 
2 red galaxies. 

(iii) At high-2, the modest difference between K-f- As and 
blue galaxies in environment distribution suggests that ei- 
ther the quenching probability is a function of environment, 
prefering overdense to underdense environments, or that 
quenching only happens to a subpopulation of blue galaxies 
which are preferentially found in overdense environments. 

(iv) At low-2, K-|-As have an environment distribution 
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indistinguishable from that of the blue galaxy population. 
This indicates that quenching probability has only a weak 
dependence on environment today. 

(v) Although we are unable to directly address the in- 
cidence of K+As in clusters due to the paucity of clusters 
in the DEEP2 volume, the existence of large K+A popu- 
lation in the field at both \ow-z and high-jz indicate that 
cluster-specific K-fA formation mechanisms cannot be the 
dominant route by which these galaxies are formed. 

Throughout this work, we have also learned a dearly- 
bought lesson of relying on measurements of K-l-A fraction 
as a function of environment. Because of the large uncertain- 
ties of environment measurements and the rarity of K-|-A 
galaxies, the significance of a K-l-A fraction vs. environment 
trend is easily overestimated and can change significantly 
when bin definitions are altered. More importantly, these 
trends are very sensitive to the composition of the all-galaxy 
sample out of which the fraction is defined. It is likely that 
these issues have led to the widely varying and sometimes 
contradictory K+A fraction-environment relations found in 
the literature. Therefore, we warn on blindly making com- 
parisons of these trends between samples defined using dif- 
ferent limiting bands, different luminosity cuts, with differ- 
ent selection rates for red and blue galaxies, or at different 
redshifts. As a result, we recommend comparing the full en- 
vironment distributions of different samples, using robust, 
nonparametric statistical tests, rather than focusing on un- 
reliable measurements of K-l-A fractions. 
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